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Abstract 

A given Hamiltonian matrix H with real spectrum is assumed tridiagonal and 
non-Hermitian, H ^ H> . Its possible Hermitizations via an amended, ad hoc 
inner-product metric Q = > are studied. Under certain reasonable as- 
sumptions, all of these metrics G = G(if) are shown obtainable as recurrent 
solutions of the hidden Hermiticity constraint G = G H called Dieudonne 
equation. In this framework even the two-parametric Jacobi-polynomial real 
and asymmetric N— site lattice H^ N \fi, u) is found exactly solvable at all N. 

1 Introduction 

Given a set 

rn(exp.) rn(exp.) rp{exp.) /-.x 

H/ , ■ • ■ , -£>at_i \l) 

of experimentally determined real eigenvalues of a quantum observable f) (one 
may think, say, about excitation energies of a heavy nucleus), we are often 
interested in its simplified theoretical interpretation. Thus, typically [lj, one 
replaces the diagonalization of the realistic (and complicated) Hamiltonian 



operator f) = f)^ by the diagonalization of its simplified isospectral "effective" 
image 

h = n-^n ^ (2) 

In the similar half-phenomenological applications, some of the most frequent 
sample sets (CQ) even carry convenient nicknames. For example, if the levels 
are equidistant, one speaks about a "vibrational spectrum". The reason is 
that we may always choose then the alternative operator H in a particularly 
simple form of harmonic oscillator. For similar reasons, the name "rotational 
band" is being assigned to the non-equidistant sets of energies Et p - ] which 
grow, with the level-numbering subscript n, quadratically. 

In the series of papers [2] - [6] we proposed an extension of the latter 
"effective-operator" philosophy to a next family of analytically tractable 
models where the values ([T]) form a less elementary series but where they 
still exhibit a certain regularity by being fitted, with a reasonable precision, 
by a friendly iV— plet 

j-i(theor.) rp(theor.) rp(theor.) /o\ 

h/ Q ,...,2j n _ x ' [6) 

defined as zeros of a suitable classical orthogonal polynomial of degree N, 



Y N (Ef neor - > ) = 0, j = 0, 1, . . . , N - 1, cf. Table dJ 

Table 1: Special tridiagonal Hamiltonians (j3J) related to classical orthogonal 
polynomials 



matrix elements (k = 1, 2, . . .) 


polynomials 


ref. 
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Formally, the exceptional character of the latter assumption resulted from 
the reality, tridiagonality and asymmetry of the N by N matrix form of the 
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underlying observable, 



a 1 ci 















b 2 a 2 c 2 











H = HW 



6 3 a 3 c 3 



(4) 
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In our present paper we intend to make the picture complete by showing how 
the construction may be made feasible for any classical orthogonal polynomial 
including the most complicated two-parametric family of Jacobi polynomials 
for which our original strategy of construction (as proposed in Ref. [2] and 
used in the other references of Table [1]) failed. 

Our present key idea will lie in the replacement of the computer-assisted 
symbolic-manipulation method of Ref. [2] by its new alternative which proves 
virtually computer- independent. In essence, we shall achieve our goal by 
making an explicit use of the tridiagonality of our classical-polynomial- related 
matrices H^ N ' and, moreover, of the positivity of their non-main diagonals. 

The presentation of our message will start, in section [2} by the discus- 
sion of a few formal aspects of the situation in which the observables of a 
given quantum system are represented by manifestly non-Hermitian finite- 
dimensional matrices. In the subsequent section [3] we shall clarify the com- 
patibility of such an approach with the standard principles of Quantum Me- 
chanics, recommending that one should rather rewrite the currently used 
but slightly misleading term "non-Hermitian" as "hiddenly Hermitian" alias 
"cryptohermitian" [7]. This hidden Hermiticity property of H is also char- 
acterized there by the so called Dieudonne equation which matches matrix 
H with the so called metric, i.e., operator = Q^Q ^ I related to Eq. fl2]). 

The key message of our present paper is delivered in section H] where the 
Dieudonne equation is shown solvable via recurrences whenever the input 
matrix H is assumed tridiagonal. The case of Jacobi polynomials is then re- 
called for illustration - this also completes the classical orthogonal polynomial 



3 



list as displayed in Table [TJ 

In the subsequent sections [5] and [6] we study, in more detail, some con- 
sequences of a pentadiagonal-matrix choice of G and H, respectively. One 
should add that while the former results illustrate the merits of the present 
recurrent approach, the latter section rather points at its natural limitations 
because, as we shall show, not all pentadiagonal Hamiltonians H ^ W can 
be made cryptohermitian via a tridiagonal metric. 

The last section [7] is summary. 

2 The doublets of Schrodinger equations 

Within the overall approach of Ref. [2] we still have to emphasize that the 
input quantities ^ and their polynomial- zero fit (j3J) need not necessarily 
represent just energies. In the full-real-line case of Ref. [5], for example, 
these values were treated as a discrete, non-equidistant grid-point quantum 
representation of the observable H = Q of a one-dimensional coordinate. 
This being said, we shall still speak here, for the sake of brevity, just about 
a "toy-model Hamiltonian H" . 

Our present identification of the secular polynomial Y^{E^ theor ^) with one 
of the classical orthogonal polynomials has a motivation in their simplicity. 
The recurrent method itself remains applicable to a more general class of 
tridiagonal-matrix models in Quantum Mechanics. In such a broadened per- 
spective, typically, one might like to fit a given sample of numerical values of 
bound state energies via any prescribed full-matrix form of H. In a prepara- 
tory step, we then merely have to employ the Lanczos tridagonalization 
method [8] and convert the given Hamiltonian into an infinite-dimensional 
tridiagonal matrix H^°°\ This matrix should further be truncated to yield its 
iV— dimensional version of the form Indeed, once our real and tridiagonal 
Hamiltonian matrix enters the linear algebraic Schrodinger equation 

H^\^ n ) = E n \^ n ), n = 0,l,...,N-l (5) 

any standard numerical technique might be used to solve it in general case. 
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The specific advantage of our present preference of having matrix (j3J) 
generated by the recurrences between classical orthogonal polynomials may 
be seen in the resulting immediate knowledge of the eigenkets \if n ) in closed 
form, 



(lVVi)) 2 



(l^n))iv-l 
V (l^n»jv J \ 



Y {E n ) 
Y x {E n ) 



\ 



Y N - 2 (E n ) 
Y N -t{E n ) J 



(6) 



At the same time, the main shortcoming of the approach may be identified 
with the related manifest asymmetry (i.e., apparently, non- Hermit icity) of 
our special matrices H^ N ' (for illustration, cf. Table [1] once more). This 
asymmetry, in particular, implies the necessity of an additional solution of 
the conjugate linear algebraic Schrodinger equation which determines the 
same real eigenvalues but the entirely different eigenvectors (which will be 
denoted, in what follows, as "ketkets" [7]), 



[HW] f |^ m )) = Em |VU , rn = 0, 1, . . . , N - 1 



(7) 



The construction of ketkets is both important (La., it completes the definition 
of a biorthogonal basis) and challenging. After all, only in the simplest case 
of Tschebyshev polynomials of paper [4J we were able to write down the 
explicit ketket analogue of the explicit kets ([SD- We shall extend the latter 
result in what follows. 



3 Two alternative interpretations of ketkets 

In a way explained in [7] and [H] the most important properties of the ket 
solutions \ip m ) of Eq. <^ and of the ketket solutions \f> m )) of Eq. (J7J should 
be seen in their mutual orthogonality alias biorthogonality 

((ip m \ipn}^0 iff m = n (8) 
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and in their biorthogonal-basis completeness, 



N-l 

Their definition is ambiguous. As long as both Eqs. (JSJ) and (J7J) are homoge- 
neous, we may choose arbitrarily normalized "input" eigenvectors \ipn) and 
\ipn )) with the real overlaps oo n = (^ipn\Wn) 7^ 0. Next, we may write down 
the general solutions of Eqs. (EJ and ((JJ) in the respective rescaled forms 

|Vfe) = |V#> X On , W = |^ ] ))x/3 n . (10) 

We will always choose the real proportionality constants a n and /?„ in such 
a manner that ((^ n |^ n ) = 1, i.e., = ^2n=o 1^™) ((^n|- This means that 
we shall only eliminate the rescaling ambiguity, say, in the kets, 

az n = az n (/3 n ) = — — . (11) 

Pn W n 

The variability of (3 n (i.e., our freedom of rescaling the ketkets) survives. 
From the point of view of its physical interpretation, it may have two forms. 
They have to be discussed separately. 

3.1 Fixed— f) scenario 

Whenever Eq. (J2J) contains a given non-unitary operator Q (so that also the 
metric becomes unique and nontrivial, = Q^Q 7^ I), the initial Hermiticity 
property of f) = f)^ proves equivalent to the easily derived cryptohermiticity 
relation 

[iz-W] + e = effW, e:=nta. (12) 

We may re- write Eq. (fl2|) as an explicit definition of the conjugate Hamilto- 
nian H>. Once we assume the reality and non- degeneracy of the spectrum 
of energies, the second Schrodinger Eq. (JTj) may be compared with the first 
one implying the proportionality rule 

|Vn}> = 7n©|V>n}- (13) 
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Constraint ((^ n \i/; n ) = 1 eliminates the apparent new freedom since 

i 1 1 2 B 2 

where function a n = a n (/3 n ) was taken from Eq. ( TTTj) . We may conclude that 
even if the metric is given in advance, the ketkets still remain ambiguous and 
may vary with the N— plet of unconstrained rescaling parameters j3 n . 



3.2 A variable— F) scenario 

Whenever we start from a given Hamiltonian and from a fixed related 
biorthogonal basis, we may invert the above procedure and reinterpret rela- 
tions (I12p as a set of constraints imposed, by the requirement of the Hermitic- 
ity of a hypothetical t), upon the metric G. Relations (fT2|) may be then called, 
for historical reasons, Dieudonne equations [TO]. They will serve us as an im- 
plicit definition determining the Hamiltonian-adapted and k— multiindexed 
family of eligible metrics = 8^ = Q^ N \H^ N \ k). Formally we obtain the 
necessary explicit formula for them by multiplying identity OH]) by matrix 
from the left. In the light of Eq. f[T3j) this yields 

N-l 

This is a finite sum which defines the general iV— parametric solution of the 
Dieudonne equation. One has to conclude that the factorizations = Q^Q 
will finally realize the correspondence between our unique input Hamiltonian 
and its multiple eligible isospectral avatars rj = f)(re). 
In this context is is worth emphasizing that each of these avatars may 
be accompanied by a few other operators g of some other observable quan- 
tities, with the corresponding (and now, evidently, multiindex-dependent!) 
pullbacks G( n \k) accompanying the original, k— independent Hamiltonian 
H^ N >. In this sense, the multiindex k carries an explicit additional physical 
information about the dynamical contents of the model (cf. a more thorough 
discussion of this point in []]). 
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Vice versa, one may weaken or even remove these multiindex-related am- 
biguities of the physical predictions via na explicit assignment of an ob- 
servable status to some of the members of the family of operators G^ n '(k). 
The best known example (and also, in parallel, one of the mathematically 
simplest ones) is the requirement of the observability of a "charge" as con- 
jectured by Bender, Boettcher and Jones [H]. In this sense, the band-matrix 
(i.e., tridiagonal, pentadiagonal etc) constructions of 9 (as used in papers 
listed in Table [1]) represent an alternative strategy which has been proposed 
and supported by several phenomenological arguments in Ref. |12j . 



4 Recurrent constructions of the band-matrix 
metrics 



4.1 Dieudonne equation and diagonal metrics 



Relation ([12]) will be treated here in the spirit of paragraph 13.21 i.e., as a 
linear set of algebraic equations, with the input given by the matrix elements 
ax, bi, . . . of our tridiagonal toy- model Hamiltonian fll]), and with the output 
giving the (non- unique [1]) definition of the matrix elements of the eligible 
metrics B = 0(H). 

As long as our Hamiltonian fll]) is a real and tridiagonal matrix, the 
simplest possible matrix form of the metric may be assumed diagonal and 
real, 

2 



e 



'N-l 







N 



(16) 



By this ansatz the Dieudonne's equation gets converted into the difference 
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Q = H^Q — H = = between the real tridiagonal matrix 





o 2 u 2 


(j 




i) 
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cA 


a 2 6 2 















c 2 6 2 


a 3 6> 3 

















C N _ 2 0N-2 


flAT_l^JV-l 


b N 0N 















a N 0N 



and its transposed partner 



o 





0iai 0ici 

# 2 &2 ^2«.2 6*2C2 

#363 9 3 a 3 9 3 c 3 "•• : 

'■■ '■• '■• '■• 

... 0N-ibN-i On-io.n-i On-iCn-i 

... OjybN Qno-n 

The diagonal elements of this difference vanish so that we are left with the 
single sequence of the recurrence relations 



(17) 



6n+ib n -\ 



n= 1,2,...,JV-1. 



(19) 



We only have to verify that all of the matrix elements in Eq. (Tl6|) remain 
strictly positive, 9j > for all j = 1, 2, . . . , iV pQ. In the diagonal-matrix 
case such a test is trivial. 



4.2 Dieudonne equation and the tridiagonal metrics 

Let us slightly simplify the conventional notation of Eq. (J3J) and redefine 



H = HW 



an a\2 

a 2 \ a 22 a 23 

a 32 a 33 a 3 i 

a 43 044 



(20) 
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postulating simply that ajk = for j > N or k > N. This will simplify the 
formulae in the first nontrivial scenario in which we replace the diagonal- 
matrix ansatz ( II 6p by its first nontrivial, tridiagonal-metric alternative 









bn 


bu 








bi2 


b 2 2 


b 2 3 








b 2 3 


h 3 




















bN-2N 














J N-2N-1 









N-2N-1 UN-1N-1 UN-1N 



'NN 



(21) 



The same recurrent approach as above may be applied when, mutatis mutan- 
dis, we only replace Eqs. ( JTT1) and ( 1181) by their appropriate pentadiagonal- 
matrix analogues. 

The main diagonal of the difference Q vanishes so that the discussion has 
to involve just the two separate diagonals of Q, viz., the outermost diagonal 
with 



Qnn+2 — a n+ln b n +ln+2 ~~ ^nn+1 a n+ln+2 ~ , Tl — 1,2, 

and the remaining condition 

Qnn+l Q"nn ^nn+l^^n+ln+l *^n+ln b nn (3 nn ^i b nn j-i dri+ln+l , 



(22) 



n 



1,2,. 



(23) 



In the first step let us turn attention to the former relations ( 1221) which define, 
in terms of an initial value, say, b\2 = 1, the sequence of the off-diagonal 
matrix elements in the metric 

, bnn+l O'n+ln+2 



Ai+ln+2 



n 



Q"n+ln 



1,2, 



N-2. 



In contrast to the brute-force symbolic-manipulation constructions as sam- 
pled in full detail in Ref. 0, the occurrence of an exceptional subscript 
n e = N such that a ne+ i ne = will not make the construction more compli- 
cated. On the contrary, we get an easier recurrence at n = N — 1. Note also 
that an artificial choice of an "illegal" 6 12 = would merely return us back 
to the previous formulae for diagonal metric. 
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In the second step let us remind the readers about the linearity of the 
Dieudonne's equation. One might make use of the knowledge of the diagonal- 
matrix solution 6 diagonal and try to work with an a— dependent tridiagonal- 
metric ansatz 

© © diagonal & Qtridiagonal ■ (p 1 ^) 

In the nondiagonal component this would enable us to select bn = without 
any loss of generality In this spirit our second, remaining set of recurrences 
( j23|) degenerates to the recurrent relations 



7 (^nnbnn+1 b nn d nn ^.\ ^n+ln+1 -in t\t i 

On+ln+1 = , n= 1,2, ...,N - 1 

^n+ln 

which define the diagonal matrix elements of Qtridiagonai- The algebraic part 
of our task is completed. One can easily verify that our recurrences reproduce 
our older tridiagonal-metric results derived by the symbolic-manipulation 
techniques and presented in the references which are listed in Table [D 

A closer inspection of the Table reveals that the symbolic-manipulation- 
extrapolation techniques were only able to cover the non-parametric cases 
(cf. the last three items 3. - 5.) or, with certain much more serious but 
still tractable technical difficulties, also the classical orthogonal polynomials 
with a single variable parameter (cf. the first two items 1. - 2.). Interested 
readers may find a more explicit description and/or a deeper analysis of a 
suppression of these difficulties in Ref. [T3] . 



4.3 The "missing example" of Jacobi polynomials 

The brute-force treatment of the "missing" Jacobi's two-parametric classical 
polynomials Y n (z) = P^' u \z) seemed to be, from the point of view of the 
extrapolation technique of Ref. [2], prohibitively complicated. At the same 
time, precisely these complications inspired our search for a new approach. 

From the point of view of applicability of our present recurrent approach, 
there is virtually no difference between the simpler and more complicated 
concrete forms of the input matrix of the (tridiagonal) Hamiltonian. One 
only has to verify the assumptions. In the case of Jacobi polynomials there 
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emerge in fact no specific problems. One just has to find the real and non- 
degenerate N— plet {-En} of roots of the N— th Jacobi polynomial, 







n 



0,1,2,. ..,N-1. 



(25) 



For any practical purposes such a secular equation may be considered "solv- 
able" because the values of E n may be found not only by means of standard 
numerical algorithms but also using their combination with many available 
specific approximation formulae [14) . 

The presence of the two auxiliary free parameters //, v > —1 may be ex- 
pected to make the resulting spectrum sufficiently flexible and, up to certain 
degree, sufficiently universal. Having in mind a spectral-fitting applicability 
of formula (I25p . it makes sense to appreciate the compact iV— dimensional 
form of the kets 



/ 



(IV>n»l 

(IV>»» 3 



\ 



(|V>n})jV-l 



\ 



\ -2 (Pn) 
(En) J 



p(M.") 
\ r N-l 



(26) 



These ket-vectors may be further compactified using the ad hoc change of 
parameters with E n — 1 := £, fi + k := fi^ and fi + v + k := a^. This yields 
the "optimal" parametrization 

/ 



1 



\ 



\ (jv-l)! ^1^2 • • • /Mf- 



+ . . . + 



CiVCAf+l • • • C"27V-2 ^ / 



(2AT-2)!! 

Obviously, it would not make too much sense to display here the explicit 
forms of the matrix elements of the metrics. These formulae would be just 
too long for the standard capacity of the printed pages of a Journal. Still, 
due to their recurrent nature, the use of standard software (e.g., of MAPLE) 
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would still keep any needed non-numerical or numerical manipulation with 
these elements feasible and straightforward. 

Having these constructions in mind (and introducing also one additional 
abbreviation v + k := u k ), we may make Table [T] complete by providing 
the explicit Jacobi-polynomial-related definition of the three diagonals of 
elements in matrix (j4]), 

a k = H k<k = , k = 1,2, . . .,N, 

&2k-20'2k 
<J k 

c k = H Kk+ i = -2k , k = 1,2, N - 1, 

b k+ i = H k+1>k = -2 ^ k , k = 1, 2, . . . , N - 1 . 

0"2fcO"2fc+l 

This makes us able to convert the recurrences for Jacobi polynomials into 
the twin Schrodinger-type eigenvalue problem where our candidate for the 
Hamiltonian is the real and tridiagonal matrix H with the desirable proper- 
ties of its outer diagonals. 

We are quite close to the climax of the construction. In particular, the 
resulting recurrences remain elementary in the diagonal-metric case where 
we obtain 

9 n+1 = 9 n cjb n+1 = 9 n kaka2k+1 , n = i, 2 , . . . , N - 1 . (27) 

Similarly, one can proceed in the tridiagonal metric case where we also just 
have to perform the insertions. For the above-mentioned typographical rea- 
sons, these insertions will be left to the readers as an exercise. 



5 Pentadiagonal metrics and the emergence 
of a cut-off-dependence 

With the pentadiagonal real and symmetric-matrix ansatz for the metric 

O ] diagonal ® ^tridiagonal ft Q pentadiagonal (28) 
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the only unknown matrix elements occur in its updated pentadiagonal part 







fell 


bu 


bis 







bi2 


b22 


^23 ^24 






bn 


b 2 3 


&33 










b 2 4 




bN-2N~l 


bN-2N 






• • bN-2N-l 


bN-lN-l 


bN-lN 







b^^2N 


bN-lN 


bNN 



(29) 



where, again, the convenient initializations 613 = 1 and 612 = &n = may 
(though need not) be used. 

The related heptadiagonal Dieudonne's equation Q = H^Q — H = is 
again just an identity along its main diagonal, Q nn — at all n — 1, 2, . . .. 
The strictly three nontrivial series of conditions have to be satisfied, therefore, 
viz., conditions 

Qnn+3 — , Qnn+2 = , Qnn+l = , U — 1,2, ... . 

The first rule 

Qnn+3 — b n+ \ n+ 3 a n +\n ~ ^nn+2 O n +2n+3 = , 71 — 1, 2, . . . 

degenerates again to the trivial recurrences for the outer diagonal of the 
metric, 

, _ b nn+ 2 d n+ 2n+3 _ -, 

On+ln+3 — , U — I, Z, . . . . 

dn+ln 

From the second rule 

bnn+2 &nn + On+ln &n+ln+2 — &nn+l On+ln+2 ~ &nn+2 0"n+2n+2 = 

we extract the recurrences which specify the intermediate diagonal of the 
metric, 

, bnn+2 0>nn ~ b nn +l &n+\n+2 — &nn+2 ®>n+2n+2 



M+ln+2 — 



0>n+lr 



, n = l,2,.... 



Finally, the third rule 

bn—ln+l Q"n—ln ~\~ &nn bnn+l ~l~ ^n+ln+1 On+lra" 
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—bnn a nn+l ~ &nn+l a n+ln+l ~ &nn+2 a n+2n+l — , 71 — 1, 2, . . . 

(where we have to add the formal definition of a i = 0) implies that we can 
readily evaluate also the remaining missing elements lying along the main 
diagonal of the metric, 

7 "n— ln+1 Q"n—ln Q"nn bnn+l 

On+ln+1 — r 

Q"n+ln 

, b nn dnn+l + b nn+ i dn+ln+1 + ^rm+2 a n+2n+l -, 

+ , 71=1,2,.... 

An insertion of the concrete input matrix elements of the Hamiltonian as 
listed in Table [T] confirms the compatibility of our older results with these 
recurrences. It is worth recalling, in particular, that the use of recurrent 
formulae reveals why the older method encountered particular difficulties via 
an emergence of the cutoff-dependence of some elements of the metric in 
multidiagonal cases. 

The pentadiagonal metrics offer the first illustration of this phenomenon. 
One may recall, e.g., Ref. [5] where an irregularity has been detected in 
the cut-off dependent pentadiagonal-metric "ultimate" element bjy^. In the 
present language, such an irregularity appears as a very natural consequence 
of the truncation. Indeed, the related recurrent definition at n + 1 = iV 
contains the component b N _ 1N+ i a N+1N which vanishes "anomalously", i.e., 
due to the truncation of the Hamiltonian cln+in = or, equivalently, of the 
metric, &at_iat + i = 0. 

Naturally, an extension of these considerations to the heptadiagonal and 
higher band-matrix metrics is straightforward and may be left to the readers. 

Another remark must be added concerning the requirement of the pos- 
itivity of the metrics. As already noticed in the preceding papers of this 
series, this is a more or less purely numerical problem, and no news emerged 
when we simplified the algebraic constructions of matrices 9. 
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6 Pentadiagonal Hamiltonians and the tridi- 
agonal metrics 

In place of a tridiagonal non-Hermitian N by N matrix Hamiltonian 
of Eq. (jlj) one might now feel inclined to extend the scope of the method 
and to try to perform the systematic recurrent reconstruction of the metrics 
= Q(H) (say, with one, three, five, . . . nonvanishing diagonals) in more 
general, multidiagonal-Hamiltonian perspective. 

Unfortunately, the same approach leads to the new, less user-friendly 
mathematical phenomena is such Let us therefore mention a few of 

them now. For the sake of definiteness let us only concentrate on the case of 
the next-to-nearest-neighbor interactions and assume that they are described 
by the following pentadiagonal real Hamiltonian 



rr(N) 
±± penta 



an 


«12 


«13 
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«23 


&24 
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«32 
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a 3 4 
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(IN-3N-1 











0>N~2N-3 


&N-2N-2 


&N-2N-1 


CIN-2N 












a N-lN-2 


a N-lN-l 















&NN-2 


OiVJV-1 


(Inn 



For the sake of brevity let us only consider the tridiagonal ansatz (121]) for 
the metric and list just a few most essential consequences of the insertion of 
these ansatzs in Dieudonne's Eq. (fT2l . 

Firstly, in a complete parallel with our preceding considerations we may 
rewrite Eq. (fT2"|) as a condition imposed upon the heptadiagonal matrix Q = 

G — G H which is real and symmetric. 

The simplest condition concerns the outermost diagonals and reads 

frn+2,ra+3 a n+2,n ~ ^n,n.+l 0'n+l,n+3 = 0, 71 = 1, 2, ... . 

This relation enables us to treat bi 2 and 623 as arbitrary initial values of 
recurrences determining all of the remaining off-diagonal matrix elements of 
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the tridiagonal metric 9, 

; b n n+1 On+l,n+3 r _ . „ 

&n+2,n+3 = , 72=1,3,5,... Or 71 = 2,4,6,.... 

In the similar manner the condition of annihilation of the next two outermost 
diagonals of Q = reads 

bn+ln+2 a n+ln + b n+ 2n+2 a n+2n ~ b nn d n n+2 — b nn +l a n +l n +2 = 

and determines the sequence of the odd or even diagonal matrix elements 

, —bn+ln+2 0, n +ln + b nn d nn +2 + b nn +l 0, n +ln+2 



J n+2n+2 



, 71 = 1,2,... 



a n+2n 

depending on the other pair of the arbitrary initial-value quantities bn and 
622, respectively. 

In the last step we are still left with the "redundant" sequence of the 
recurrences 

bn+ln+2 O-n+ln+l + b n+ 2n+2 &n+2n+\ + b n +2n+3 a n +3n+l — 
— b nn+ i a nn +2 — bn+ln+l a n+ln+2 ~ ^n+ln+2 a n+2n+2 = , TL = 0, 1, 2, ... . 

In contrast to the tridiagonal-Hamiltonian results we are now merely allowed 
to set 6-1,0 — at n = 0. Obviously, we are left with the set of the formally 
superfluous but still perfectly valid constraints. Thus, these conditions must 
be read as a restriction imposed directly upon the input Hamiltonian itself. 
In this role these recurrences read 

^n+2rt+2 d n +2n+l + b n +2n+3 a n+3n+l 

0>n+2n+2 ~ O-n+ln+l H 



bn+ln+2 

bnn+1 dnn+2 + &n+ln+l ^n+ln+2 n , 

7 , 71 = 0,1,2,... 

O n +ln+2 

and act, with the fifth arbitrary initial value an, as a mere step by step 
recurrent definition of the diagonal matrix elements of the "admissible" or 
"tri diagonally hermitizable" pentadiagonal input Hamiltonian H^ n \ a . 
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7 Conclusions 



We can summarize that it is now possible to declare the classical-polynomial 
oscillators exactly solvable. In our present paper, our recent proposal [2] 
of the simulation of a given empirical set of energy levels (or, alternatively, 
of an N— plet of discrete experimental eigenvalues of any other quantum 
observable) by the N— plet of zeros of a suitable classical orthogonal poly- 
nomial Yn(E) found its final formulation based on the vector-structured, 
recurrent treatment of the matrix Dieudonne equation. The menu of pos- 
sible sufficiently elementary identifications of an iV— plet of components of 
wave functions with the polynomial sequences Y n (E), n = 0, 1, . . . , N — 1 has 
been made complete. 

Originally, our present study has been motivated by the failure of the 
assignment of the metrics to the "last missing" Jacobi-polynomial oscillators 
by the extrapolation method of Ref. [2]. A new approach has been devel- 
oped, in which the three-term recurrence relations which are satisfied by the 
classical orthogonal polynomials were redesigned as finding their main use in 
Dieudonne equation itself. 

Although the existence of the three-term recurrence relations satisfied by 
the classical orthogonal polynomials Y n (x) is a very well known fact, the use 
of these polynomials in the role of components of a wave function always 
seemed to be hindered by the manifestly non-Hermitian appearance of the 
Schrodinger-equation resembling truncated versions of these recurrences. 

The situation started to be changing when people realized that the ver- 
ification of Hermiticity is in fact strongly dependent on our selection of the 
inner product of wave functions. A true breakthrough (and a slow acceptance 
of the new paradigm) emerged when the manifestly non-Hermitian appear- 
ance of the imaginary cubic potential has been recognized as "reparable" (cf., 
e.g., the details of the history as outlined in [15]). 

In such a context the appeal of the classical-orthogonal-polynomial-related 
tridiagonal real and asymmetric Hamiltonians appeared to be twofold. 
Firstly, their finite-dimensional "effective-matrix" form fitted very well the 
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mathematical requirement of having an exactly solvable quantum model, 
without too many formal subtleties and with a nontrivial metric given in 
closed form. Secondly, the physical information happened to be carried in 
these models, in a not too usual manner, by the doublet of operators, viz, by 
a given H and deduced 0. 

The double appeal of our H^s in the context of model-building might 
find its further support in potential applications in the domain of quan- 
tum lattices and chain models. Two minor obstacles are encountered and 
removed. Firstly, in the context of mathematics one gets rid of the appar- 
ent non- Hermit icity of the natural tridiagonal candidate H ^ W for the 
Hamiltonian via an ad hoc, non-unitary (a.k.a. Dyson's) isospectral map Q. 
Secondly, in the context of physics, an additional merit is that the set of 
the eligible metrics is restricted to the band-matrix ones, characterized by a 
partial preservation of the concept of the localization of the lattice sites [3J. 

The key technical advantage of the model has been found in the tridiago- 
nality of its Hamiltonians. The attempted move beyond the family of classical 
polynomials (i.e., the attempted transition to the pentadiagonal and further 
Hs) has been found to lead to a loss of simplicity. For the tridiagonal Hamil- 
tonians, on the contrary, we have got rid of the apparent non-Hermiticity 
of the most natural tridiagonal candidates H ^ W for the Hamiltonians in 
full generality. This was achieved via the use of the concept of a non-unitary 
(a.k.a. Dyson's) isospectral map Q connecting a certain unknown operator 
\) (which is assumed Hermitian but, presumably, prohibitively complicated) 
with our friendly and explicitly diagonalizable matrix H^ N \ Moreover, the 
set of the eligible metrics = Q^Q has been successfully ordered as starting 
from a short-ranged sub-hierarchy characterized by the (2k + 1)— diagonal 
matrix structure of with k <C [(N + l)/2]. Thus, at the small k — 0, 1, . . ., 
a partial return to the concept of a "smeared" locality of the lattice has been 
achieved [3j y] . 

In our present paper we stressed the importance of making the assign- 
ment of a menu of metrics to a preselected Hamiltonian H less dependent on 
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the simplicity of the elements of H. In all of our older papers on the subject 
(cf. their list in Table [T|) such a simplicity has been required as a necessary 
condition of the applicability of the truncation-extrapolation method as pro- 
posed (and successfully illustrated) in Ref. [2J. With the new method we 
were now able to complete the list and incorporate also the "last missing" 
two-parametric Jacobi polynomials into the overall scheme. 

The reasons of the success of the construction may be found in a cer- 
tain recurrent structure which has been found hidden in the N by N matrix 
form of the Dieudonne equation. Unfortunately, the one-directional recurrent 
nature of the construction H — > Q{H) is only preserved in the tridiagonal- 
Hamiltonian cases. The situation becomes less favorable for the pentadi- 
agonal and higher generalizations of H since one then must fine-tune the 
compatibility of the metric G with Hamiltonian if in a selfconsistent man- 
ner. 
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